Quantum control of the hyperfine-coupled electron and nuclear spins in alkali atoms 
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We study quantum control of the full hyperfine manifold in the ground-electronic state of alkali 
atoms based on applied radio frequency and microwave fields. Such interactions should allow essen- 
tially decoherence-free dynamics and the application of techniques for robust control developed for 
NMR spectroscopy. We establish the conditions under which the system is controllable in the sense 
that one can generate an arbitrary unitary on the system. We apply this to the case of ^^^Cs with 
its d = 16 dimensional Hilbert space of magnetic sublevels in the 6S1/2 state, and design control 
waveforms that generate an arbitrary target state from an initial fiducial state. We develop a gener- 
alized Wigner function representation for this space consisting of the direct sum of two irreducible 
representations of SU{2), allowing us to visualize these states. The performance of different control 
scenarios is evaluated based on the ability to generate a high-fidelity operation in an allotted time 
with the available resources. We find good operating points commensurate with modest laboratory 
requirements. 

PACS numbers: 32.80. Qk,42.50.-p,02.30.Yy 



I. INTRODUCTION 

The control of spins is the foundation of coherent spec- 
troscopy at the heart of NMR, atomic clocks, and many 
precision metrology experiments [ ] . More recently, spins 
have been seen as ideal carriers of information, with de- 
velopments in spintronics for classical [ ] and quantum 
[3] information processing. Atomic spin systems have 
been of particular interest given their excellent isolation 
from the environment and the available techniques in 
the "quantum optics toolbox". Examples include en- 
sembles of atomic spins as quantum information process- 
ing elements [4, 5, 6, 7, 8], ion-trap quantum computers 
[9, 10, 11], and neutral-atom optical lattices [ ]. The lat- 
ter has attracted tremendous attention in recent years, as 
controllable spin lattices are seen as a platform in which 
to perform quantum simulations of condensed matter sys- 
tems [LI] and studies of topological quantum field theory 
[14]. 

While in many studies of atomic spin control one con- 
siders two- level spin qubits, real atoms have large spins 
with a rich internal structure. The ability to fully control 
the Hilbert space within the atoms for various applica- 
tions is an important addition to the toolbox. It allows 
for the possibility of d-dimensional qudits as the funda- 
mental information carriers [ ] and the embedding of 
logical qubits in a qudit, which may be advantageous for 
control or protection from errors [ ]. Additionally, ma- 
nipulating a nontrivial Hilbert space allows us to explore 
interesting dynamics such as quantum chaos [ , ] . Fi- 
nally, in the same way that liquid state NMR has pro- 
vided an excellent platform for exploring quantum con- 
trol protocols [19, 20, 21], atomic spin systems provide a 
test-bed with unique physical properties that allows for 
new investigations into control and measurement tech- 
niques. 

In this paper we study quantum control of electron and 



nuclear spins of alkali atoms, coupled by the hyperfine in- 
teraction in the electronic ground state, using combina- 
tions of static, AC-radio-frequency, and AC-microwave- 
frequency magnetic fields. In previous studies, we imple- 
mented similar control based on a combination of mag- 
netic interactions and a nonlinear AC-Stark shift induced 
by a laser field [22, 23]. In that work, control was re- 
stricted to a single subspace of total angular momentum 
of the coupled spin system, rather than the whole Hilbert 
space. More fundamentally, the light-shift interaction at 
the heart of the protocol came at the cost of some deco- 
herence by spontaneous emission. The maximum ratio of 
nonlinear light shift to photon scattering is fixed by the 
atomic structure, thereby limiting the ultimate utility 
of that approach. In contrast, direct magnetic coupling 
to spins in the ground state is essentially decoherence 
free, with dephasing due solely to inhomogeneities and 
background fields that can be mitigated, in principle, by 
robust control techniques [2(J, 24]. Rf/ microwave con- 
trol thus has the potential for higher-fidelity operation 
on a larger Hilbert space with speeds comparable to or 
faster than those previously achieved. Such capabilities 
are similarly being explored for use in ion trap quantum 
information processing [25]. 

Our ultimate goal is the implementation of general dy- 
namical maps on the quantum system. In this paper we 
take a first step - preparation of an arbitrary state in the 
Hilbert space. In particular, we look at open-loop state 
preparation through the application of control waveforms 
that take some particular known fiducial state to a target 
final state. It was shown in [' ] that this type of problem 
is much easier than the problem of generating arbitrary 
unitary evolutions, due to a mathematical property that 
promises local searches will find global optima. We will 
discuss this in more detail in Sec. IV, but intuitively this 
problem requires specifying only one row of a d x (i uni- 
tary matrix. Creating one row requires as many free pa- 



rameters as it takes to identify a single state in Hilbert 
space, 2(i — 2, as opposed to the cP — 1 parameters that 
are required to specify an arbitrary unitary matrix. 

The remainder of this article is organized as follows. 
In Sec. II we establish the fundamental Hamiltonian that 
describes the dynamics of the spins and their interaction 
with the applied control fields. In Sec. Ill, we discuss 
what it means to be controllable in the Lie algebraic sense 
[27, 28, 29] and determine which configurations of our 
fields satisfy these criteria. Finally, in Sec. IV we present 
an algorithm that uses a gradient ascent method to gen- 
erate control waveforms for state preparations. We apply 
our protocol to the example of ^^^Cs, with one valence 
electron (spin S = 1/2) and a nuclear spin (/ = 7/2), 
a total Hilbert space of dimension d = 16. We show 
data from simulations comparing how different types of 
controls perform and tradeoffs that can be expected in 
the laboratory. Along the way, we establish new visual- 
ization tools for coupled spin systems based on a gener- 
alized Wigner function representation and give analytic 
proofs for controllability of our system. These results are 
discussed in the appendices. 



II. CONTROL HAMILTONIAN 

In this work we seek to control the quantum state of 
a multilevel atom. Though single atom addressing and 
measurement are possible [30, 31, 32], in practice we con- 
sider ensembles of uncorrelated particles. To the degree 
that the atoms are identically prepared and uniformly ad- 
dressed, with no interactions between them either from 
interatomic forces or through measurement backaction, 
we can take the joint state of the system as effectively N 
identical copies, p^^ . More general many-body control 
is not considered here. Restricting then to a single atom, 
the relevant Hilbert space of an alkali atom in its elec- 
tronic ground state is the tensor product space of elec- 
tronic spin S and nuclear spin / subsystems, H = i)s^^i- 
Given the single valence electron S = 1/2, the Hilbert 
space is spanned by two irreducible subspaces of total an- 
gular momentum F± = / ± 1/2, such that 7i = f)+ ()_. 

The Hamiltonian describing the atom and its interac- 
tion with external magnetic fields takes the form. 



H = Al-S^2pBB{t)-S, 



(1) 



where ps is the Bohr magneton and we have neglected 
the small nuclear magneton contribution. Here and 
throughout we set fi = 1. We consider the application 
of three fields, B{t) = Boe^^Brf{t)^B^^{t). The static 
bias field Bq defines the quantization axis and Zeeman 
splittings between the magnetic sublevels. The terms 
Brf(t) and B^^{t) describe magnetic fields oscillating at 
radio and microwave frequency, respectively. The hyper- 
fine coupling between spins provides an effective nonlin- 
earity that will allow full controllability of the Hilbert 
space for appropriate choice of external fields. 



In the linear Zeeman regime, /hbBq <C A^ the static 
field acts separately in the two irreducible subspaces, and 
according to the Lande projection theorem the Hamilto- 
nian is approximately. 



Hbo ^ MB ^ g/Bo • F^^) 



(2) 



/=± 



Here F^^^ = P±FP± refers to the total angular momen- 
tum operator projected onto the subspaces with quantum 
number F±. Neglecting the nuclear magneton contribu- 
tion, the g-factors for the two manifolds have equal mag- 
nitude but opposite sign, i.e. ^+ = —g- = 1/F+. The 
hyperfine coupling plus bias magnetic field thus deter- 
mine the static Hamiltonian, 



-no - -7, (•f + 



P_) + l^o(Fi+)-i^i-)), (3) 



where AEhf = AF-^ is the hyperfine splitting and Qq = 
jibBq/Fj^ is the Zeeman splitting between neighboring 
magnetic sublevels. 

As our first control field, we consider rf-magnetic fields 
oscillating near the frequency of the Zeeman splitting, 
cjrf ~ r^o, realized by Helmholtz coils driven with the 
appropriate current. We take two sets of coils that pro- 
duces fields with x and y polarization, independent am- 
plitude and phase control, but equal carrier frequency, 
a;rf. Again, for moderate current such that the ampli- 
tude of the field is in the linear Zeeman regime, the rf- 
Hamiltonian takes a form equivalent to the interaction 
with the static field 

H,i{t) = n,{t) cos {ojrit - Mt)) {Fi+'> - Fj-)) 

+ ny{t) cos {urit - <Py{t)) (F^+) - F^-)). (4) 



F = 4 



F = 3 




FIG. 1: The ground state hyperfine manifold of ^^^Cs. Rf- 
magnetic fields (in red) lead to independent rotations in the 
two manifolds. Microwaves (in blue) are the generators of 
rotation in a two-dimensional subspace between states in the 
two manifolds, here the stretched state transition |4,4) -^ 
|3,3). 



The time dependent amplitudes ((^^^(t), (l^yC^)) ^^^ 
phases {(l)x{t),(j)y{t)) of the two sets of rf coils will be 
used to control the system. 

To better understand the effect of the rf field, con- 
sider a resonant interaction, cjrf = ^^o- ^^ the ro- 
tating frame, Hj.^{t) -^ H^^{t) = Uj^HYf{t)Urfj where 
[/rf = exp < —iUrft{Fz — Fz ) MS a rotation of the two 

manifolds about the z-axis in opposite directions, Fx -^ 
Fi^^ cos(cjrft) ± F^^^ sin(cjrft), F^^^ -^ F^^^ cos(cjrft) T 
Fx sin(c(;rf^). Performing this unitary transformation 
and averaging over a cycle, the rf-Hamiltonian in the ro- 
tating wave approximation is. 



H'At) 



^cos(0.(t))(Fi+)-Fi- 



sin(</..(i))(FW+F(-)) 
M)eos(0,(t))(FW-F(-)) 
^sin(^,(t))(Fi+)+Fi-)) 



(5) 



Rf-control of the two spin manifolds differs from the fa- 
miliar spin resonance problem. In the latter, a single 
magnetic field in either the x or ^/-direction would be 
sufficient to generate the entire SU{2) algebra for rota- 
tions. With two irreducible manifolds there is an added 
freedom - the two angular momenta F+ and F_ can ro- 
tate in the same or opposite directions. Amplitude and 
phase control of two rf-magnetic field polarizations al- 
lows us to perform arbitrary and independent rotations 
on the two hyperfine manifolds. With only a single direc- 
tion of Brf we would be restricted to either co-rotating 
or counter-rotating in the two subspaces. 

The weak rf-magnetic fields alone will not be sufficient 
to fully control our atomic system; they don't couple the 
F+ and F+ manifolds, nor do they provide a nonlinear 
Hamiltonian within these subspaces. In order to make 
our system fully controllable, we look to resonant mi- 
crowaves. While the fundamental Hamiltonian governing 
the microwaves is exactly of the same form as the qua- 
sistatic magnetic fields, the resonant behavior leads to 
very different dynamics than the previous interactions. 
Depending on the polarization and frequency, the mi- 
crowave couples a Zeeman sublevel in F+ manifold with 
one in the F_ manifold whose magnetic quantum num- 
ber differs by Am±0, 1. For a sufficiently strong bias Bq 
we can ignore any off-resonant excitation, and restrict 
the Hamiltonian to act only on a 2D subspace spanned 
by the states we are trying to couple. In that case the 
microwave Hamiltonian has the form 



H^^{t) = ^^w{t) cos (cj^w^ - (l)^w{t))ax, 



(6) 



where ax is the Pauli sigma-x matrix for this pseudospin, 
ax = |F+,m+)(F_,m_| + |F_,m_)(F+,m+| andn^^{t) 
is the (time-dependent) Rabi frequency depending on 
the microwave power and the transition matrix element. 



Again, the amplitude and phase of the microwave fields 
are control parameters. In this subspace, the problem 
takes the form of the standard two-level resonance prob- 
lem. We must take care in going to the rotating frame to 
account for the simultaneous transformation we perform 
due to the rf-fields. The complete frame transformation 
is achieved by the unitary 



U = /7rf exp 



at , _ 



P-) 



(7) 



where a = cj^^ — (^+ + m_)cc;rf. Under this transfroma- 
tion, the Hamiltonian in the rotating wave approximation 
for resonant microwaves is 



H'{t) = %^cos(0,^(i))<T. 



^h^sm{4>,^{t))a 



y 



(8) 



generating rotations of this pseudo-spin on the Bloch 
sphere. 

Combining the static, rf, and microwave interactions 
the final Hamiltonian in the rotating frame is 



H'{i) 



% + H[,{t) + H'^^{t). 



(9) 



Allowing for a finite detuning of the oscillating fields from 
resonance, the static Hamiltonian in the rotating frame 
becomes. 



H', = ^{P, 



P_) + A,f(Fi+)-Fi-)), 



(10) 



where A 



/iW 



■^/iW 



AEhf — {m-\- m')ix)^i is the effective 
detuning of the microwaves from the two-level transition 
of interest, |F_,m_) -^ |F+,m+), and Arf = uo^i — ^o is 
the rf detuning. This, together with Eqs. (5,8), defines 
the Hamiltonian we employ for control, and which we 
will analyze for use in arbitrary state preparation. 



III. CONTROLLABILITY 

In order to perform state preparation on a system, we 
must first determine the conditions under which the sys- 
tem is controllable in principle. Is it possible to create an 
arbitrary state using the Hamiltonian dynamics we have 
available, neglecting technical constraints such as band- 
width and slew rates that restrict the types of waveforms 
we use to drive our atoms? To answer this, it is simpler to 
analyze the conditions necessary to generate an arbitrary 
unitary evolution, a problem that has been studied in 
depth in the control theory literature, [27, 28], and more 
recently from a quantum information perspective [2y] . 

Formally stated, we consider a quantum system in a 
Hilbert space of dimension d^ governed by a Hamiltonian 
of the form 



H{t) = Ho + Y,h,{t)Hj. 



(11) 
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FIG. 2: Table exploring controllability of the system for a variety of configurations: one microwave field driven on different 
two-level transitions, \F — 3, M) ^ \F — 4, M^), amplitude and/or phase control, one or two sets of orthogonal rf coils (rf 
polarizations), and resonant vs. detuned fields. The different configurations yield one of four different outcomes: (green circle) 
all microwave transitions provide full controllability, (yellow square) all transitions but the clock transition |3,0) -^ |4,0) 
provide full controllability, (orange pentagon) only the transitions of the form |3,±3) -^ |4, ±4) and |3,±3) -^ |4, ±2) provide 
full controllability, and (red octagon) no transitions yield controllable Hamiltonian dynamics. In this calculation we consider 
all valid microwave transitions that can be selected with polarization and/or frequency. 



The system is said to be "controllable" , if for every pos- 
sible unitary map, Uq G SU{d)^ there exists a choice of 
controls bj (t) and a finite time T such that the Hamilto- 
nian evolution given by the Schrodinger equation, U{t) = 
—iH(t)U{t)^ maps the identity operator to Uq at time 
T. A necessary and sufficient condition for controlla- 
bility is that the independent terms in the Hamiltonian 
{Hq^ i^i, . . . , Hn} generate the Lie algebra Qu{d). Obvi- 
ously, if we can perform any unitary evolution, we can 
perform any state preparation. 

To elucidate the connection between controllability 
and the generators of the Lie algebra, we review the ba- 
sic principles here. A Lie algebra is a linear vector space 
with an algebraic product defined by the commutator. 
We can see that we can generate any linear combination 
of our initial set of generators by looking at very short 
square-pulses according to the Trotter formula, where 



tions it is also possible to generate the commutators by 



-iHiaA^-iH2f3A 



-iiaHi+f3H2)A 



(12) 



Such short pulses are allowed since we assume access to 
an arbitrary waveform. In addition to linear combina- 



^-iHiA^-iH2A^iHiA^iH2A 



-[Hi,H2]A^ 



(13) 



The ability to generate, in principle, any linear combi- 
nation and any commutator means that one can simu- 
late any element of the the Lie algebra generated by our 
initial independent Hamiltonians, {Hq^ il^i, . . . , Hn}^ and 
thus any unitary contained in the associated Lie group. 
If the Lie algebra generated is 5u((i), we call this system 
controllable. 

For the Hamiltonian system described in Sec. II, 
with arbitrary control of the amplitude and phase of the 
two orthogonal sets of rf-coils and a single microwave 
field, the control algebra generated by the six operators 
{F^'^\F^^\F^~\F^~\(Ta,, dy} is 5u{d) in its entirety. In 
this case, it is possible to prove controllability analyt- 
ically for an arbitrary alkali, with an arbitrary nuclear 
spin /. The proof is fairly involved and is shown in detail 
in Appendix A. 

Though sufficient, the entire available set is not neces- 
sary to achieve controllability. In practice, one can reduce 



the number of generators in the control algebra and still 
implement an arbitrary unitary. For an experiment, it 
is important to understand which components are really 
necessary so that we can evaluate the tradeoffs between 
ease of implementation and controllability. In order to 
study the capability of various reduced sets of controls 
we resort to numerics. Being a linear vector space, deter- 
mining whether a set of operators generates the algebra 
5u{d) only requires showing that it is possible to generate 
a basis for 5u{d). We do this by tabulating a library of 
all the linearly independent operators generated by the 
set of control operators. 

We carried out this procedure for the specific example 
of ■'^^^Cs with nuclear spin / = 7/2 and study the capa- 
bility of a variety of control sets to generate the entire 
5u(16) algebra. We considered 8 different microwave con- 
figurations: controlling or fixing the amplitude and the 
phase of the fields, and whether or not we are detuned 
from resonance. The two cases where both the ampli- 
tude and phase are controlled and where the amplitude 
is fixed but the phase is controlled can be shown to be 
equivalent. In the rf configurations, we also allow for one 
or two orthogonal sets of magnetic coils (rf polarization). 
The last free parameter is the choice of which microwave 
transition we excite. We assume arbitrary frequency and 
polarization selectivity of the desired transition for this 
purpose. The results are summarized in Table 2. In 
each box we enumerate the set of microwave transitions 
that yield controllable dynamics. We find that our sys- 
tem is controllable for a wide number of configurations, 
though there are some specific cases in which it is not. 
For example, out of all the choices for microwave tran- 
sitions, the clock transition, |i^+,0) -^ |F_,0), renders 
the system controllable in the least number of scenarios. 
This should not come as much of a surprise since we are 
controlling the system with rf magnetic fields and this 
transition is insensitive to magnetic fields. 

It is interesting to note that there exist configurations 
that are controllable in which there is one time-dependent 
control field and some fixed time-independent fields. This 
is the simplest scenario one could expect to find, and 
allows for bang-bang control, a well-studied protocol. In 
this paper, however, we look at the control systems that 
utilize more parameters, decreasing the time needed for 
state preparation. 



IV. STATE PREPARATION 



We seek to design Hamiltonian evolutions that take 
an initial known quantum state to an arbitrary quan- 
tum state in the Hilbert space. Historically, numerical 
searches for control waveforms have performed much bet- 
ter than expected. The fidelity of a state preparation is 



a functional of the control waveform, 

F[h{t)] = \{^target\U[h{t)]\4>o)\' 

(14) 

where T is the time ordering operator. Under ideal 
conditions, assuming no decoherence and an arbitrarily 
amount of time to perform the control, Rabitz et al. [26] 
proved that the control landscape is surprisingly simple 
- every local optimum is a global optimum. This implies 
that, 

Vb(.)i^[bo(t)] = ^ F[bo(t)] = {0, 1}, (15) 

completely independent of the initial and target states. 
Therefore, a local search of the space of control fields, 
starting from any random initial guess, will find a global 
maximum of the fidelity. For this problem, gradient 
searches perform about as well as more computationally 
intensive searches like genetic or simulated annealing al- 
gorithms. 

In a real system the assumptions of the proof will not 
hold. There will always be some decoherence and one 
does not have infinite time to perform the control. In 
fact, we would like to perform state preparation as fast 
as possible in order to combat decoherence and various 
inhomogeneities that lead to accumulated errors. Addi- 
tionally, we need to consider control fields that have a 
limited bandwidth and slew rate constraints. For these 
realistic conditions, not every gradient search from an ar- 
bitrary starting point yields a global maxima. Nonethe- 
less, we have found empirically that the results of the 
theorem are approximately true with moderate decoher- 
ence and after a sufficient time. We still find excellent 
protocols after making only a small handful of searches, 
and these can be further filtered to find control waveforms 
that perform well under realistic operating conditions. 



A. Optimization protocol 

As we are dealing with the optimization of waveforms 
that are functions of continuous time, the first step is to 
transform the problem into a search for a finite number 
of values at discrete times. The physical constraints of 
bandwidths and slew rates of the controllers provide a 
natural scale. There is a minimum interval during which 
a field can vary over a maximum range. A discretized ver- 
sion of a control waveform is thus specified as a vector 
of values within this range at these fixed intervals. The 
continuous control waveforms are then found by interpo- 
lation using cubic splines, consistent with the bandwidth 
constraints, at least on a fine enough grid for use in our 
numerical integration of the Schrodinger equation. 

We create optimal control waveforms by first fixing 
the total time of the state preparation procedure. Due 
to our discretization technique, fixing the total time fixes 
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FIG. 3: Sample evolutions that result from our state preparation algorithm. Starting in the spin coherent state |4,4), we 
simulate preparation of the state 4=(|4,4) + |3, —3)) and obtain a fidelity 0.993. We control the amplitudes and phases of 
rf coils in both the x and y directions, as well as the amplitude and phase of a resonant microwave that couples the states 
|4, —4) and |3, —3). We show the Cartesian components of the three control fields (QcoscI) and Qsincf)) over the entire state 
preparation time of 150/is. We show snapshots of the evolved state at five different times, identified as times (0)-(4). Two 
different representations of the state are shown: bar charts of the absolute values of the density matrix elements, and a 
generalized spherical Wigner function. The spheres on the diagonal represent the Wigner functions in the irreducible subspaces 
F^ and the off-diagonal spheres represent the coherences between the manifolds. For details see Appendix B. 
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FIG. 4: Same as Fig. (3), preparing the state ^|4, 4) + |(|3, 3) + |3, -3)) with a fidehty of 0.995. 



the number of optimization variables. Starting from a 
randomly chosen initial vector of control waveform val- 
ues, bo, we perform a gradient ascent search by taking 
small steps in the direction of steepest ascent, i.e. 



bn+i=bn + eVF(b^). 



(16) 



An optimal value corresponds to the maximum, where 
the gradient approaches zero. We performed this search 
numerically on a Matlab cluster by optimizing waveforms 
from a handful of random seeds in parallel, and then 
chose the one that gave the highest fidelity. In this work 
we do not consider the robustness of the waveform to in- 
homogeneities and noise. More complex objective func- 
tions can be optimized as in [ ] once the relevant exper- 
imental conditions are known. 

We applied this protocol to the specific case of ^^^Cs, 
with ground-state hyperfine splitting of AEhf = 9-2 
GHz. We take a static bias field to produce a Zeeman 
splitting of r^o = 1-0 MHz, sufficient to give excellent 
resolution of the magnetic sublevels, but well within the 
linear Zeeman regime. The rf field power is chosen so 
that on resonance the rotation rate is characterized by 
r^rf = 15 kHz. As a generic case, we take one mi- 
crowave field, resonant on one of the stretched transitions 
|F = 3, M = ± - 3) ^ |F = 4, M = dz - 4), where the 
microwave Rabi frequency is largest, and the system is 
controllable in a wide variety of scenarios. The microwave 
power is chosen to give a Rabi frequency l^^w = 40 kHz. 
The slew rates constrain the maximum rate of change of 
amplitude and phase of the control fields. In the case of 
the rf-magnetic field, a "slew time" of Trf = 10/is fixes the 
slew rates on the amplitude to 1.5 kHz//is and phase to 
0.2 7r//is. In the case of microwaves, faster control is pos- 
sible, with a slew time of r^w = 1.0 /is, or amplitude and 
phase slew rates of 40 kHz//is and 2.0 7r//is respectively. 

Two examples of the end product of this optimiza- 
tion are shown in Figs. (3,4) for target states 4=(|4, 4) + 

|3,-3)) and ^|4,4) + ^(|3,3) + |3,-3)) respectively. 
The initial state for these examples is the stretched state 
1 4, 4), a state easily reached by optical pumping. We op- 
timized waveforms for both the amplitude and phase for 
two rf polarizations in the x and y directions as well as 
for the microwave field. The fidelities of preparation in 
both cases are greater than 99%. The state preparations 
shown here take 150/is, an interval that ensures that a 
moderate search will yield high-fidelity waveforms. More 
intensive optimizations can yield faster control wave- 
forms. 

Our gradient search algorithm leads to waveforms that 
cause the system to undergo quite complex dynamics, as 
evidenced by the intermediate states seen in the course 
of the evolutions. Figs. (3,4). One may wonder whether 
there are simpler choices, since given a fixed initial state, 
there are many different waveforms that lead to same tar- 
get state. While our method does lead to waveforms that 
are hard to intuitively understand, some recent studies 
[33] suggest that the waveforms derived from gradient 



searches may be more robust than those that come from 
more geometric algorithms. 



B. Performance of Optimization 

In Sec. HI we discussed the mathematical conditions 
necessary for our Hamiltonian dynamics to be control- 
lable. These conditions, while useful for ruling out large 
classes of Hamiltonians as unsuitable for our purposes, 
tell us nothing about the relative performance of differ- 
ent control scenarios. Our figure of merit is the time after 
which we can be reasonably sure that our optimization 
will find a high fidelity waveform for any target state. To 
determine this time for a given control protocol, we run 
our optimization up to a given final control time over a 
large collection of randomly chosen states and determine 
the average fidelity. In this section we examine these re- 
sults and discuss some of the tradeoffs and bottlenecks 
that might be encountered in the lab. 

There are many parameters in this system that we 
can manipulate, including the number of independently 
controlled rf polarizations, the number of resonant mi- 
crowave frequencies, the types of controls (amplitude vs. 
phase), detuning, slew rates, and the strengths of the 
different fields. Based on some of our previous exper- 
iments we set as a baseline: one microwave frequency, 
two orthogonal rf polarizations, rf power giving l^rf =15 
kHz, a microwave Rabi frequency of Q^^ = 40 kHz, a 
rf slew time of 10 /is, and a microwave slew time of 1.0 
/is. While we could independently vary all these param- 
eters, this would be an unwieldy computation. Here we 
fix some of the parameters that are unlikely to differ in 
the future experiments we are considering. In particular, 
we fix the rf slew time to be 10 /is and consider control 
with two sets of rf coils. For simplicity we also consider 
all fields to be resonant, and the microwaves to couple 
the stretched states. 

Statistics were collected by running the state prepara- 
tion algorithm for 10 different random states found by 
sampling using the Harr measure on SU{16) [)4]. In all 
cases the initial state was the |4, 4) state. For each com- 
bination of total time, target state, and system config- 
uration, we run the optimization 20 times starting from 
different random seeds of the vector that defines the con- 
trol waveform, as discussed is Sec. HI. Out of this set of 
20, we choose the highest fidelity preparation. The fideli- 
ties from the 10 random states are averaged to produce 
the data points shown in Fig. 5. In principle, more it- 
erations would yield higher fidelity waveforms, but it is 
useful to understand which types of high-fidelity controls 
can be found after only modest searches. 

In Fig. 5a, we study the effect of varying the charac- 
teristics of the microwave field. We compare the perfor- 
mance of one vs. two resonant microwave frequencies on 
one or both of the stretched transitions, |3,3) -^ |4,4) 
and 1 3, —3) -^ |4, —4). In addition, we examine the effect 
of removing control of the microwave amplitude (a see- 
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FIG. 5: Plots of the average fidelity of state preparation for different control configurations and total preparation times. Each 
point represents the fidelity averaged over a set of 10 states randomly chosen from the Harr measure. For each state and 
configuration, the gradient search was performed with 20 random seeds and we chose the protocol that generated the highest 
fidelity. 



nario that still allows for full controllability of the system, 
as discussed in Sec. III). As expected, since the microwave 
Rabi frequency is larger than the rf Larmor frequency, 
increasing the number of microwave fields has a large ef- 
fect. On the other hand, it was surprising that fixing 
the microwave amplitude, thereby substantially decreas- 
ing the number of control parameters, yielded higher fi- 
delity waveforms. We suspect that while there most likely 
exist higher fidelity waveforms with control of both am- 
plitude and phase, increasing the number of microwave 
control parameters rapidly increases the dimension of the 
search space, requiring many more iterations of our al- 
gorithm to find a superior waveforms, on average. This 
suspicion is reinforced by Fig. 5b, where we consider the 
effect of microwave slew time. With our baseline parame- 
ters, it would appear that increasing the microwave slew 
time doesn't really limit the optimized control perfor- 
mance. In fact, the smallest slew time we considered, 1.0 
/is, performed slightly worse than the other slew times, 
including the case where the microwave amplitudes are 
fixed. As a reminder, the slew time determines the infor- 
mation content of our waveforms, and thus the number of 
optimization variables. Again, we see that for the mod- 
est searches we are performing, decreasing the dimension 
of the search space counterbalances the loss of control. 

In Fig. 5c, we study the effect of the power in the rf and 
microwave fields. For these simulations we fixed the am- 
plitudes of the fields and solely control their phases. We 
find that varying the microwave power around our base- 
line makes little difference. The rf power is slightly more 
important, but increasing the Larmor frequency above 
the baseline has a fairly small effect. These results in- 
dicate that the slew rate and bandwidth constraints we 
have imposed on the rf magnetic fields are the bottleneck 
for controlling the system, and limit the ability to more 
rapidly control the system through increases in power. 



It would appear that the microwave parameters we em- 
ploy as our baseline are also well above the limits im- 
posed by this bottleneck and we can safely reduce the 
microwave power and slew rates without sacrificing per- 
formance. The rf Larmor frequency we employ is com- 
mensurate with the slew rate constraint. 

By optimizing many state preparations for a variety 
of control configurations we find state preparation proto- 
cols with this system that take between 50 — 150/is. We 
can compare this to the types of control waveforms that 
were implemented in our previous work that employed a 
nonlinear AC-Stark shift to achieve controllability [23]. 
The waveforms we find here are about an order of mag- 
nitude faster, control a Hilbert space that is double the 
dimension, and have negligible decoherence as compared 
to the intrinsic decoherence that arises from spontaneous 
emission. 



V. SUMMARY AND OUTLOOK 

In this paper we have studied quantum control of the 
d = 2(2/ -h 1) dimensional Hilbert space associated with 
coupled electron spin 5 = 1/2 and nuclear spin / of alkali 
atoms in their electronic ground state, based on interac- 
tions with static, rf, and microwave magnetic fields. Such 
interactions allow rapid and essentially decoherence-free 
dynamics. We studied a variety of configurations that 
allow for full controllability of the system based on ana- 
lytic proofs for the most general control fields considered 
and numerical studies in more restricted configurations. 
With controllability in hand, we studied the problem of 
open-loop state preparation, mapping a known fiducial 
state to an arbitrary target state of the Hilbert space, 
applied to the specific problem of ^^^Cs with a d = 16 
ground-electronic subspace. Control waveforms can be 
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found from simple gradient searches in the control land- 
scape. We evaluated the performance of a variety of sce- 
narios, restricting some control parameters by e.g. fixing 
the amplitudes of the fields or the number of resonant mi- 
crowaves frequencies. We find that under certain condi- 
tions, restricted control yielded better performance. We 
attribute this to the complexity in searching a large di- 
mensional control parameter space. 

Implementation of the proposed control protocol dis- 
cussed in this article will require diagnostics to measure 
the fidelity of the prepared state with respect to the tar- 
get. This can be done via quantum state tomography 
on the ensemble. In prior work, we developed and im- 
plemented a protocol whereby the quantum state is esti- 
mated via continuous measurement on a single ensemble 
of identically prepared atoms [*^^, ''^^]. To achieve this, 
the system must be controllable. By applying a well- 
chosen waveform in the course of the continuous measure- 
ment, one gains access to an informationally complete 
record. One then inverts the measurement history to ob- 
tain a high-fidelity estimate of the initial state, limited 
only by the signal-to-noise ratio and decoherence that 
occurs during the measurement. In our prior work, com- 
binations of laser interactions and magnetic fields were 
employed to yield an informationally complete measure- 
ment record. In that case, the laser beam acted both as 
a probe of the atoms to provide the measurement record 
and as a control field to provide a nonlinear light-shift on 
the atoms. To extend this protocol to the case at hand, 
we must control the system to produce an information- 
ally complete set of observables on the full d-dimensional 
Hilbert space. Using the microwave/rf interactions stud- 
ied here, one can achieve this while separating the control 
and measurement functions of the applied fields. This 
should be faster and reduce decoherence induced by pho- 
ton scattering of the probe. We will study this in future 
work with the goal of rapid and robust quantum state 
estimation in a large dimensional Hilbert space. 

In addition to state preparation. Lie-algebraic control- 
lability implies that there exist waveforms to generate 
arbitrary unitary maps on the system. The work con- 
sidered here corresponds to implementing one column 
of a unitarity because the waveforms we design lead to 
the intended dynamics only on a single fiducial quantum 
state. Nonetheless, such capabilities provide a starting 
point to more general control tasks such as embedding 
a qubit in a qudit, or the implementation of universal 
qudit control. Moreover, because the Hamiltonians con- 
sidered here correspond to generators of rotations either 
within irreducible subspaces or on pseudospins, it gives 
us a natural starting point to consider the application of 
composite pulse sequences developed for NMR [2u] in or- 
der to make our protocols more robust to the inevitable 
imperfections in our system. 

This research was supported by NSF Grants No. 
PHY-0653599 and No. PHY-0653631, ONR Grant No. 
N00014-05-1-420, and lARPA Grant No. DAAD19-13- 
R-0011. 



APPENDIX A: ANALYTIC CONTROLLABILITY 

We prove here that our quantum system is control- 
lable given accesses to the Hamiltonian presented in 
Sec. II. More precisely we show that the set of operators 

{Fx \Fy \Fx iFy iCFxiCTy} generates the Lie algebra 
5u((i), where d is the dimension of the tensor product 
space, d = 2(2/ + 1). In the language of control theory, 
a set of operators is said to "simulate" another operator 
if we can construct this operator through linear combi- 
nation and commutators. Starting with a generating set, 
one simulates new operators which are added to our li- 
brary. The goal is to use the generators to simulate a 
basis for the entire Lie algebra. 

We begin by proving a very general theorem relating 
to controllability when one is able to perform SU(2) ro- 
tations on an n-dimensional Hilbert space. 

Theorem 1 In an d-dimensional Hilbert space with d > 
2, if one has access to the irreducible generators of rota- 
tions, Jx and Jy, then in order to fully control the space 
it is sufficient to add an operator h that has a non-zero 
overlap (according to the trace inner product) with at least 
one rank-2 irreducible spherical tensor. That is 

3 q s.t Tr {hT^^^^ ^0 ^ { J„ Jy, h}L.A. = 5u(d). 

Here we have introduced the orthonormal basis of ir- 
reducible spherical tensor operators. 



rf H-^) = ^J^^J2^J,m+q\k,q;J,m)\J,m+q){J,m\, 

m 

(Al) 



satisfying the fundamental commutation rules. 



(A2) 



where J± = Jx :^iJy. It follows from these commutators 

that given the set {Jx-,Jy,Tq ^} one can simulate any 
rank-/c irreducible tensor, and since these are an opera- 
tor basis, the generators of rotation can map any rank-A; 
operator to any other rank-/c operator. With this prop- 
erty we are now prepared to prove a lemma. 

(2) 

Lemma 1 {Jx^Jy^T^ ^} generates 5u{d) . 
We prove this by first noting that 



7^(2) rp{k) 

^0 '^g 



Ck. 



T?+')+4aT?-'^ 



q-^q 



(A3) 



The exact form of the constants is irrelevant except for 
the fact that there is always some rank-/c tensor for which 
Ck^q is nonzero. Given this, the proof follows by induc- 
tion. Suppose our library of simulatable operators con- 
tains all operators of ranks k and k — 1. By commuting 
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some rank k operator with Tq ^ we obtain an operator 
with support on operators of rank k — 1 and /c + 1, thus 
containing a component in the space of rank k -\- 1 op- 
erators that is hnearly independent from the current set 
of Hamiltonians in our hbrary. Commutation with the 
generators of rotation ahow us to simulate ah other rank 
k -\- 1 operators. Since we can simulate all rank-1 from 
the generators {Jx^Jy}^ and the rank-0 operator is the 
trivial identity operator, it follows by induction that we 
can simulate all rank-/c operators that are supported on 
the Hilbert space, k < d — 1. Therefore {JxjJy^T^ } 
generates 5u{d). QED 

With this lemma, we see that in order to show theorem 
1, we need merely to show that the set {Jx^Jy^h} can 

(2) 

simulate the operator Tq . We will do this in essentially 
three steps. Before we start we expand the Hamiltonian 
h in our spherical basis, h = X]/c=i ^q=-k ^Q ^Q • 



1. Step 1: Simulate hi 



rp{2) , ■^-vd— 1 ■^-vfc 






To simulate hi we note that h is defined to have some 
nonzero rank-2 component. With rotations we can trans- 

(2) 

form the rank-2 component to Tq \ Additionally, since 
we have all the rank-1 tensors in our library already, we 
can remove the rank-1 piece of h through linear combi- 
nations to yield hi. 



2. Step 2: Simulate /i2 = T^''^ + J^tZl K^^^T^^^ 
Consider the double commutator 

[Jz, [Jz, h^]] =J2i2 9'<'^r('=). (A4) 

k=3q=-k 

If we take a linear combination hi — a[Jz, [Jz^hi]] the 
resulting operator has the same coefficients for g' = 0. For 
qo 7^ 0, choosing a = l/g^, we can sequentially remove 
all rank-2 tensor components, and we are left with /i2- 



3. Step 3: Simulate T^'^^ 

Consider the double commutator 



[Jx^ [Jx^ ^2], 



7^(2) 


+ -^'(Ti^^^ 


+ r!^2^) 


d-1 
k=3 


h'^'^ (2kik 


+ l)T^''^'> 



We repeat the process in Step 2 to remove the compo- 
nents from /i2 with (7 7^ to obtain 



h'o 



3^(2) 

2^« 



d-l 
k=3 



k{k- 



^0 • 



(A6) 



If we now take the linear combination /12 — 2/^2/(^0(^0 + 

1)) we remove the Tq °^ component, but are left with a 

(2) 
nonzero Tq ^ term. Repeating this procedure for ko = 

3 ... (d — 1) yields an operator that is proportional to 



n(2) 



This completes our proof of theorem 1. 



With theorem 1 in hand we are now prepared to show 
that {Fx , Fy , Fx , Fy'' .(Jx.cry} generates the Lie al- 
gebra 5u{d) . Unfortunately, we do not start off with the 
ability to perform the irreducible generators of rotations 
on the entire space, and so cannot immediately prove 
controllability using theorem 1. We do, however, have 
access to the generators of angular momentum for both 
the T+ and T_ subspaces, and so we begin by showing 
that we can simulate any operator that has support on 
only one of the two manifolds. 

To show controllability of the T+ manifold we re- 
quire an operator that has a nonzero overlap with a 
rank-2 tensor on that space. Restricted to the T+ sub- 
space, the (Tz operator looks like a projector onto some 
particular sublevel, |T+, 777,+ )(T+, 7ti+|. The overlap of 

this projector with Tq^^^ is Tr (^1^+, m+)(T+, mj,\T^^A = 

Y^5/11(T+, m+|2, 0; T+, m+), which is nonzero for all val- 
ues of 771+ . Of course, a^ has support in the T_ man- 



ifold. However, 



r(+) 



is confined to the T+ mani- 
fold. Since commuting by Fx can't change the rank of 
a tensor, we are left with an operator confined to the T+ 
manifold that has a nonzero overlap with some rank-2 
tensor, and so according to theorem 1, we have complete 
control of the T+ manifold. This proof directly carries 
over to the T_ manifold. 

At this point we have shown that we have full control- 
lability over both the T+ and the T_ subspaces, as well 
as the 2-dimensional subspace coupled by the resonant 
microwaves. We can write this in matrix form 



Sl 






Cr's 








S2 



(A7) 



where we have ordered the basis vectors so that the 
-f a/(A: — l)k(k + l)(k + 2)(T2^^ -^T^)) ). states coupled by the microwaves are adjacent to each 

/ other. We have shown that we can simulate any opera- 
(A5) tor that only has matrix elements within the three boxes 
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FIG. 6: Representations of states with bar charts of the absolute values of the density matrix elements, (a.i) is a spin 
coherent state \il^)ai = |4,4) and (a.ii) is a superposition two oppositely oriented spin coherent states, one for each of the two 
manifolds, \ip)aii — 4|(|4,4) + |3, —3). In (b.i, b.ii) we show the effects of rotations on a superposition of spin squeezed states, 

each determined as the ground state of F^ — Fy in the respective irreducible manifold. Finally, in (c.i) we have a coherent 
superposition of the state |4, 0) and a cat state ^(|3,3) + |3, —3)) and in (c.ii) we have an incoherent mixture of those two 
states. 



in Eq. A7. This includes all diagonal operators on the 
whole matrix as well as all operators with only nonzero 
matrix elements on the super-diagonals (one off the di- 
agonal). The irreducible representations of angular mo- 
mentum, Jx and Jy^ on the entire space have support 
only on the super-diagonals. Therefore, we can simulate 
Jx and Jy. According to theorem 1 all we need to show 



for controllability is that we can simulate some operator 
with a nonzero overlap with some rank 2 operator. Since 
we can simulate any diagonal operator, we can simulate 
T^^^ . It thus follows that {fJ^^ , F^'^^ , fJ"^ , F^~^ ^cTx.cTy} 
generates 5u{d). 
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FIG. 7: Representations of the six states shown in Fig. (6) by the generahzed spherical Wigner functions. Each state is 
represented by four spheres. The spheres on the diagonal are the standard SU(2) Wigner functions in the F = 4 (upper 
diagonal) and F = 3 (lower diagonal) irreducible subspaces. The radius of these spheres, ranging from zero to one, determines 
the total population in that subspace. The off-diagonal spheres represent the coherences between the two subspaces (see text). 



APPENDIX B: GENERALIZED WIGNER 
FUNCTION REPRESENTATION 



In dealing with high dimensional spin systems, it is 
useful to be able generate graphical representations of 
the quantum states which give some geometric intuition. 
The spin coherent state Wigner function representation 
introduced by Agarwal [ ] provides a generalization of 



the standard Wigner function based on harmonic oscilla- 
tor coherent states used to describe infinite dimensional 
systems. Given a spin J, the spin coherent state Wigner 
function is essentially a multipole representation on the 
sphere defined as, 



W^{e,cj^) = 5]5]Tr[pfW(J)]rf)(0,</.), (Bl) 
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where Yq {O^cj)) are the spherical harmonics, and Tq ^ (J) 
are the irreducible spherical tensors given in Eq. (Al). 
For a given spin, the indices describing non-trivial irre- 
ducible tensors run from < /c < 2J and —k<q<k. 
These plots are useful visualization tools because they 
capture the effect of geometric rotations on the quantum 
state. Two quantum states that differ solely by a SU(2) 
rotation will generate Wigner functions that also differ 
from each other by the same physical rotation. 

We seek to generalize this to the case of a tensor prod- 
uct space of two spins (here electron and nuclear), equiva- 
lent to the direct sum of two irreducible representations of 
SU(2) in the hyperfine subspaces, F and F^ We achieve 
this by considering the expanded set of tensors defined 

by 



TJ.''\F,F') 



2k + 1 



2F- 



-J2{F,m + q\k,q;F',m}\F,m + q){F',m\. 



(B2) 



The range of the indices is now \F — F'\ < k < 
F -\- F' and —k < q < k. One can easily show 
that for two spin manifolds, the set of operators 
{tP {F, F), Tf ) (F, F'), Tf ) (F', F), T^ ^ (F', F')} com- 
prises a complete orthonormal operator basis for the ten- 
sor product space. We again can map these operators 
to the spherical harmonics, and for each state get four 
spherical Wigner functions: one each for the F and F^ 
manifolds, and two for the coherences between manifolds. 
We label them Wf,FjWf,f'jWf',f and Wf',f'' By the 
Hemiticity of the density operator, Wf,f' and Wf',f con- 
tain redundant information and are complex, so one need 
only consider the real and imaginary part of Wf,f' , yield- 
ing four real functions. 

We scale the radii of the spheres over which the Wigner 
function is plotted. For the functions that describe a 
given hyperfine manifold, we let the radius of the sphere 
equal the population in the subspace, Tt{PfpPf)- In or- 



der to set the radii of the spheres corresponding to the 
coherences between the manifolds, we look at the sum of 
the singular values of the off-block component of the den- 
sity matrix, y/^rn^m' \{F^m\p\F' ,m^)\'^. This allows 
for nonequal dimensions of the two subspace. Addition- 
ally, we scale these "coherence spheres" by the ratio of 
the real versus imaginary parts of Wigner function. The 
primary purpose of doing this is to be able to distinguish 
between pure superpositions and incoherent mixtures be- 
tween the two manifolds. 

To gain some intuition, we show examples of dif- 
ferent states and different representations. Figure (6) 
shows bar charts of the absolute values of the density 
matrix elements for the six states: \^|^)a^ = |4,4) and 
\i^)aii = (|4, 4) + |3, — 3))V^ are spin coherent states and 
their superposition; \^|^)b^^ and \ilj)bii are superpositions 
of spin squeezed states in the two manifolds along dif- 
ferent quadratures; \ip)ci^ and |V^)cii are coherent su- 
perpositions vs. incoherent mixtures of a "cat state" 
(|3,3) + |3, — 3))/a/2 in one manifold, and a Dicke state 
|4,0) in the other. The corresponding Wigner functions 
are shown in Fig. (7). From these plots we make the 
following observations. When restricted to a subspace 
corresponding to a given hyperfine manifold, the Wigner 
functions on the diagonal have the familiar forms of 
SU{2) Wigner functions, with the radius of the sphere 
determining the total population in that subspace. The 
off-diagonal Wigner functions show the effect of the co- 
herences, had the entire Hilbert space been determined 
by an irreducible representation. This is clearly seen in 
Fig. (7aii), where the coherences are of the familiar form 
for a superposition of "north" and "south" pole spin co- 
herent states. The effect of geometric rotation is exhib- 
ited in |'0)5i and \tlj)bii- The bar charts do not indicate 
any similarity between the states, while the Wigner func- 
tions are clearly related by a 90 degree rotation. Finally, 
the difference between coherent superpositions and inco- 
herent mixtures of states in the two manifolds is clearly 
seen in Fig. (7c). 
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